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We present details from a recent test of Lorentz and CPT symmetry using hydrogen masers 
We have placed a new limit on Lorentz and CPT violation of the proton in terms of a recent 
standard model extension by placing a bound on sidereal variation of the F = 1, Attif = ±1 
Zeeman frequency in hydrogen. Here, the theoretical standard model extension is reviewed. The 
operating principles of the maser and the double resonance technique used to measure the Zeeman 
frequency are discussed. The characterization of systematic effects is described, and the method of 
data analysis is presented. We compare our result to other recent experiments, and discuss potential 
steps to improve our measurement. 



I. INTRODUCTION 

A theoretical framework has recently been developed that incorporates Lorentz and CPT symmetry violation 
into the standard model and quantifies their effects. [p|-^14[. One branch of this framework emphasizes low energy, 
experimental searches for symmetry violating effects in atomic energy levels In particular, Lorentz and CPT 

violation in hydrogen has been examined and sidereal variations in the = 1, Ami? = ±1 Zeeman frequency have 
been quantified Motivated by this work, we have conducted a search for sidereal variation in the hydrogen 

Zeeman frequency, and have placed a new clean bound of 10~^^ GeV on Lorentz and CPT violation of the proton [Q. 
Here we provide additional details of the theoretical framework, experiment and analysis. In Sec. || we discuss the 



standard model extension. In Sec. HI we describe the basic concepts of hydrogen maser operation and our Zeeman 
frequency measurement technique. In Sec. |^ we describe the procedure used to collect data and extract a sidereal 
bound on the Zeeman frequency. In Sec. ^ we describe efforts to reduce and characterize systematic effects. Finally, 
in Sec. Vl we compare our result to other clock-comparison tests of Lorentz and CPT symmetry, and discuss potential 



means of improving our measurement. 

II. LORENTZ AND CPT SYMMETRY VIOLATION IN THE STANDARD MODEL 

Experimental investigations of Lorentz symmetry provide important tests of the standard model of particle physics 
and general relativity. While the standard model successfully describes particle phenomenology, it is believed to be the 
low energy limit of a fundamental theory that incorporates gravity. This underlying theory may be Lorentz invariant, 
yet contain spontaneous symmetry-breaking that could result in small violations of Lorentz invariance and CPT at 
the level of the standard model. 

A theoretical framework has been developed to describe Lorentz and CPT violation at the level of the standard 
model by Kostelecky and coworkers This standard-model extension is quite general: it emerges as the low- 

energy limit of any underlying theory that generates the standard model and contains spontaneous Lorentz symmetry 
violation |^-^. For example, such characteristics might emerge from string theory A key feature of the standard 

model extension is that it is formulated at the level of the known elementary particles, and thus enables quantitative 
comparison of a wide array of searches for Lorentz and CPT violation ||9[-p^ . 

"Clock comparison experiments" are searches for temporal variations in atomic energy levels. According to the 
standard model extension considered here, Lorentz and CPT violation may produce shifts in certain atomic levels, 
whose magnitude depends on the orientation of the atom's quantization axis relative to a fixed inertial frame ■ 
Certain atomic transition frequencies, therefore, may exhibit sinusoidal variation as the earth rotates on its axis. 
New limits can be placed on Lorentz and CPT violation by bounding sidereal variation of these atomic transition 
frequencies. 

Specifically, the description of Lorentz and CPT violation is included in the relativistic Lagrange density of the 
constituent particles of the atom. For example, the modified electron Lagrangian becomes |l3] 



1 



where 



and 



£ = iz^r.a> - W + cZ'' (1) 



M = m + l^a^Y + ^^757'^ + ^^m-'^''') ■ (3) 

The parameters a^, 6^, c^^, c?^^, e^, /i^, g\^i, and iJ^;^ represent possible vacuum expectation values of Lorentz 
tensors generated through spontaneous Lorentz symmetry breaking in an underlying theory. These are absent in the 
standard model. The parameters a^, e^, and ^a/ji^ represent coupling strengths for terms that violate both 
CPT and Lorentz symmetry, while Cfn,, dfj_^, and H^j_i, violate Lorentz symmetry only. An analogous expression exists 
for the modified proton and neutron Lagrangians (a superscript will be appended to difi^erentiate between the sets 
of parameters). The standard model extension treats only the free particle properties of the constituent particles, 
estimating that all interaction effects will be of higher order As a result, the interaction term Cf^^^ is unchanged 
from the conventional, Lorentz invariant, QED interaction term. 

Within this phcnomcnological framework, the values of these parameters are not calculable; instead, values must 
be determined experimentally. The general nature of this theory ensures that different experimental searches may 
place bounds on different combinations of Lorentz and CPT violating terms, while direct comparisons between these 



experiments are possible (see Table ID and Ref. |13|). 

The leading-order Lorentz and CPT violating energy level shifts for a given atom are obtained by summing over the 
individual free particle shifts of the atomic constituents. From the symmetry violating correction to the relativistic 
Lagrangian, a non-relativistic correction Hamiltonian Sh is found using standard field theory techniques p3[ | . Assuming 
Lorentz and CPT violating effects to be small, the energy level shifts are calculated perturbatively by taking the 
expectation value of the correction Hamiltonian with respect to the unperturbed atomic states, leading to a shift in 
an atomic [F^mp) sublevel given by 

AEF,mf. = {F,mF\ne5he + UpShp + nnShn\F,mF)- (4) 

Here riu, is the number of each type of particle and Sh^ is the corresponding correction Hamiltonian. Note that 
for most atoms, the interpretation of energy level shifts in terms of this standard model extension is reliant on the 
particular model used to describe the atomic nucleus (e.g., the Schmidt model). One key advantage of a study in 
hydrogen is the simplicity of the nuclear structure (a single proton), with its results uncompromised by any nuclear 
model uncertainty. 

Among the most recent clock comparison experiments are Penning trap tests by Dehmelt and co-workers with the 
electron and positron |p^ , |l7| which place a limit on electron Lorentz and CPT violation at 10^^^ GcV. A recent 
re-analysis by Adclbcrgcr, Gundlach, Heckel, and co-workers of existing data from the "Eot-Wash 11" spin-polarized 
torsion pendulum has improved this to a level of 10^^^ GeV the most stringent bound to date on Lorentz 

and CPT violation of the electron. A new limit on neutron Lorentz and CPT violation has been placed at 10"'^^ GeV 
by Bear et al. using a dual species noble gas maser and comparing Zeeman frequencies of ^^^Xe and '^He. The 
current limit on Lorentz and CPT violation of the proton is 10^^^ GeV, as derived from an experiment by Lamoreaux 
and Hunter which compared Zeeman frequencies of ^^^Hg and ^■^^Cs. 

Figure |l| shows the Lorentz and CPT violating corrections to the energy levels of the ground state of hydrogen . 
The shift in the F = 1, Amp = ±1 Zeeman frequency is p3[: 



\Ai^z\ = l\iK- dl,m, - H^u) + {b^3 " ^fo^p - Hf^)]. (5) 



The subscripts denote the projection of the tensor couplings onto the laboratory frame. Therefore, as the earth 
rotates relative to a fixed inertial frame, the Zeeman frequency vz will exhibit a sidereal variation. We have recently 
published the result of a search for this variation of the -F = 1, Ato^^ = ±1 Zeeman frequency in hydrogen using 
hydrogen masers . This search has placed a new, clean bound on Lorentz and CPT violation of the proton at a 
level of 10-27 GeV. 
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FIG. 1. Hydrogen hyporfino structure. The full curves are the unperturbed hyperfinc levels, while the dashed curves illus- 
trate the shifts due to Lorentz and CPT violating effects with the exaggerated values of I63 — d^o'm,,. — Hl2\ = 90 MHz and 
I63 — df^^rrip — Hf2\ = W MHz. This work reports a bound of less than 1 mHz for these terms. A hydrogen maser oscillates on 
the first-order magnetic field-independent |2) ^ |4) hyporfino transition near 1420 MHz. The maser typically operates with a 
static field less than 1 mG. For these low field strengths, the two F = 1, Amp = ±1 Zeeman frequencies are nearly degenerate, 
and z/12 « 1/23 ^ 1 kHz. 
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III. HYDROGEN MASER CONCEPTS 



The electronic ground state in hydrogen is spht into four levels by the hyperfinc interaction, labeled (following 
the notation of Andresen [Q) |1) to |4} in order of decreasing energy (Fig. |l]). The energies of atoms in |1) and 
|2) decrease as the magnetic field decreases; these are therefore low-field seeking states. Conversely, |3) and |4) are 
high-field seeking states. In low fields, |2) and |4) are only dependent on magnetic field in second order. The maser 
oscillates on the |2) ^ |4) transition (field-independent to first-order). This transition frequency, as a function of 
static field, is given by 1^24 = i^hfs + 27505^ {v in Hz with B in Gauss, with v^fs ~ 1420.405751 MHz the zero- 
field hyperfinc frequency). Hydrogen masers typically operate with low static fields (less than 1 mG), where V24 is 
shifted from Vfifg by about 3 mHz, or 2 parts in 10^^. The two F — 1, Amp — ±1 Zeeman frequencies are given 
by V12 = 1.4 X lO'^B - 13755^ and V23 = 1-4 x lO^B + 13755^ At B = 1 mG these are nearly degenerate, with 
1^12 — '-'23 ~ 3 mHz, much less than the Zeeman linewidth of approximately 1 Hz. 



A. Maser operation 

In a hydrogen maser P5|-p7[ , molecular hydrogen is dissociated in an rf discharge and a beam of hydrogen atoms 
is formed, as shown in Fig. A hexapole state selecting magnet focuses the low-field-seeking hyperfinc states |1) 
and |2) into a quartz maser bulb at about 10^^ atoms/sec. Inside the bulb (volume ~ 10^ cm"^), the atoms travel 
ballistically for about 1 second before escaping, making ^ 10"* collisions with the bulb wall. A Teflon coating reduces 
the atom-wall interaction and thus inhibits decoherence of the masing atomic ensemble by wall collisions. The maser 
bulb is centered inside a cylindrical TEqh microwave cavity resonant with the 1420 MHz hyperfinc transition. The 
microwave field stimulates a small, coherent magnetization in the atomic ensemble, and this magnetization acts as 
a source to stimulate the microwave field. With sufficiently high atomic fiux and low cavity losses, this feedback 
induces active maser oscillation. The maser signal is inductively coupled out of the microwave cavity and amplified 
with an external receiver. Surrounding the cavity, a solenoid produces the weak static magnetic field (~ 1 mG) that 
establishes the quantization axis inside the maser bulb and sets the Zeeman frequency (« 1 kHz). A pair of Hclmholtz 
coils produces the oscillating transverse magnetic field that drives the F = 1, Amp = ±1 Zeeman transitions. The 
cavity, solenoid, and Zeeman coils are all enclosed within several layers of high permeability magnetic shielding. 

A well engineered hydrogen maser can have fractional stabilities approaching 10~^^ over intervals of hours. This 
stability is enabled by a long atom- field interaction time (1 s), a low atom- wall interaction (due to the low atomic 
polarizability of H and the wall's Teflon coating), reduced Doppler effects (the atoms are confined to a region of 
uniform microwave field phase, effectively averaging their velocity to zero over the interaction time with the field), 
and multiple layers of thermal control of the cavity (stabilizing cavity pulling shifts). 



B. Maser characterization 



Among the quantities used to characterize a hydrogen maser, those most relevant to this experiment are the atomic 
line-Q Qz, the population decay rate 71, the hyperfinc decoherence rate 72, the atomic flow rate into and out of the 
bulb 7b, and the maser Rabi frequency IX24I. We describe here a comprehensive set of measurements to characterize 
hydrogen maser P-8. The results discussed here are summarized in Table |. Our Lorentz and CPT symmetry test 
data were taken with a similar but newer hydrogen maser, P-28 ^8|] . A few of the maser characterization parameters 
for P-28, while not directly measured, have been inferred using fitting parameters from the double resonance method 



used to measure the F = 1, Amp — ±1 Zeeman frequency, described in Sec. |IV A| . These values are included in Table 
I in italics. 

To determine these parameters of an operating H maser, the cavity volume Vc, bulb volume Vf,, cavity quality 
factor Qc, filling factor 77, and output coupling coefficient (3 must be known. For both masers, Vc = 1.4 x 10~^ m'^, 
Vb = 2.9 X 10-3 m3, Qc « 40,000, and /3 = 0.23 ||. The filling factor, defined as |§ 

^ = -igilM^, (6) 



cavity 



quantifies the ratio of average magnetic field energy inside the bulb to the average magnetic field energy in the cavity. 
This has a value of 77 = 2.14 for masers P-8 and P-28 [pOl. 
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FIG. 2. Hydrogen maser schematic. The solenoid generates a weak static magnetic field Bo which defines a quantization axis 
inside the maser bulb. The microwave cavity field Hp (dashed field lines) and the coherent magnetization M of the atomic 
ensemble form the coupled actively oscillating system. 



5 



parameter 


symbol 


P-8 


cavity volume 


Vc 


1.4 X 10"^ m^ 


bulb volume 


Vt 


2.9 X 10"=* m^ 


cavity-Q 


Qc 


39,346 


filling factor 




2.14 


line-Q 


Qi 


1.6 X 10^ 


maser quality parameter 


q 


0.100 


maser relaxation rate 


7t 


1.83 rad/s 


bulb escape rate 


lb 


0.86 rad/s 


population decay rate 


71 


4.04 rad/s 


maser decoherence rate 


72 


2.77 rad/s 


spin-exchanee decay rate 


7s e 


1.06 rad/s 


radiated power 


P 


600 fW 


threshold power 


Pc 


250 fW 


output coupling 


/9 


0.23 


output power 


Po 


112 fW 


total flux 


Itot 


15.0 X 10^^ atoms/s 


flux of |2) atoms 


I 


3.13 X 10^^ atoms/s 


threshold flux 


Ith 


0.54 X 10^^ atoms/s 


atomic density 


n 


2.8 X 10^^ atoms/m^ 


maser Rabi frequency 


\X2A 


2.77 rad/s 



P-28 



2.14 
1.6 X lO'' 



0.86 rad/s 
2.88 rad/s 
2.8 rad/s 



0.23 
100 fW 



2.14 rad/s 



TABLE I. Maser chara cteriz ation parameters. The italicized values for P-28 were inferred from double resonance flt param- 

AU other values were either calculated or extracted from direct measurements as described in 



eters as described in Sec. fVA 
this section. 



For canonical hydrogen maser operation, there are two important relaxation rates . For a room temperature 

H maser, the decay of the population inversion is described by the longitudinal relaxation rate 

7l =76+7r + 27se+7l, (7) 

and the decay of the atomic coherence is described by the transverse relaxation rate 

72 = 7b + 7r +7se +72- (8) 

Here, 7b is the atomic flow rate into the bulb, 7^ is the rate of recombination into molecular hydrogen at the bulb wall, 
7se is the hydrogen-hydrogen spin-exchange decay rate, and 7- includes all other sources of decay, such as decoherence 
during wall collisions and effects of magnetic field gradients. 

In the steady state, the atom flow rate into the bulb is equal to the geometric escape rate from the bulb, given by 
7h — vA/AKVb, where v = 2.5 xlO^ cm/s is the mean thermal velocity of atoms in the bulb, A ~ 0.254 cm"^ is the 
area of the bulb entrance aperture, and K Ki 6 is the Klausing factor [Q. Thus, jh — 0.86 rad/s for both P-8 and 
P-28. The spin exchange decay rate is given approximately by |^ , p9| 

1 _ 

7se = -nvrcr (9) 

where Vr — 3.6 xlO^ cm/s is the mean relative velocity of atoms in the bulb and cr = 21 x 10^^^ cm^ is the 
hydrogen-hydrogen spin-exchange cross section. The hydrogen density is given by p^j29| ] 

(10) 



(lb +lr)Vb 

where Itot is the total flux of hydrogen atoms into the storage bulb. 

The atomic line-Q is related to the transverse relaxation rate and the maser oscillation frequency w by ||2£ 

272 



(11) 
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It is measured using the cavity pulling of the maser frequency: neglecting spin-exchange shifts, the maser frequency 
is given by pH 



iO = UJ24 + — (i^C - ^24) • 



(12) 



By measuring the maser frequency as a function of cavity frequency setting, the line-Q can be determined. For both 
P-8 and P-28, we find Q; = 1.6 x 10^, and therefore 72 = 2.8 rad/s. 

A convenient single measure of spin-exchange-independent relaxation in a hydrogen maser is given by "gamma-t" 
P7|,p9| 



7* = [(76 + 7r + 7^(76 + 7r + 72) 



(13) 



Using this, a more useful form for the longitudinal relaxation rate, 71, can be found. By combining Eqn. |l3| with 
Eqns. and ||, we find 



71 



7t 



■ 27se- 



72 -7s 

Using Eqns. BO, we can relate the line-Q to /, the input flux of atoms in state |2) as |2£ 



1 



, 1 

7fc + 7r + 72 + 9i— 7t 



Qi 

using the threshold flux required for maser oscillation (neglecting spin-exchange) 



Ith — 



47r /ilQc??' 



and the maser quality parameter 



aVrh 



It 



76 +7r 



Vc_ 
vVb 



1 



hot 



(14) 



(15) 



(16) 



(17) 



The ratio I /hot is a measure of the effectiveness of the state selection of atoms entering the bulb. While / is not 
directly measurable, it can be related to the power P radiated by the atoms by P7|,P9| 



P 



+ (1 - 3<z) 



where Pc = htolth/'^- The maser power is also related to the maser Rabi frequency by E? 



P = 



inuj 1x24 1 2 



7172 



IX 



24 I 



7172 



(18) 



(19) 



The power coupled out of the maser is given by Po/P — P/i^ + P)- 

Generally, the parameter q is less than 0.1, while I/Ith is approximately 2 or 3. Hence, the first term of Eqn. |l^ 
can be neglected relative to the others. If we make the reasonable approximation that 7^ =72, then we can rewrite 
Eqn. |l5| using Eqns. p|, |l6Hl| as (29| 



= mP + b 



(20) 



where 



, 2 

b = —It 



(21) 



and 



7 



1 

It 



(22) 



Therefore, by measuring the hne-Q as a function of maser power and extracting the slope m and the y-intercept 6, 
we can determine q and 7f. For maser P-8, q = 0.100 and 7f = 1.83 rad/s. 

With these values of q and 7t, we found Ith — 0.54 x 10^^ atoms/s (using Eqn. |6|), and Pc = 250 fW. With a 
measured output power of Po = 112 fW, the atoms were radiating P — 599 fW, and the flux of state |2) atoms was 
/ = 3.13 X 10^^ atoms/s (Eqn. |l^). Under the assumption that 7t ~ 76 + 7r we found that the total flux was Itot = 
15.0 X 10"'^^ atoms/s (Eqn. |lj) and the density was n = 2.8 X 10^5 atoms/m^ (Eqn. |l|). The spin-exchange decay 
rate was then found to be ^se = 1-06 rad/s (Eqn. Finally, the population decay rate was 71 = 4.04 rad/s (Eqn. 
and the maser Rabi frequency was \X2i\ = 2.77 rad/s (Eqn. |l9|). 



C. Zeeman frequency determination 

The F = I, Amp = ±1 Zeeman frequency is measured using a double resonance technique p^ , |3l| , ^ . As the 
frequency of an audio frequency magnetic field w^, applied perpendicular to the quantization axis, is swept through 
the Zeeman frequency, a shift in the maser frequency is observed (Fig. ^) . When the applied field is near the Zeeman 
frequency, two-photon transitions (one audio photon plus one microwave photon) link states |1) and |3) to state |4), 
in addition to the single microwave photon transition between states |2) and |4). This two photon coupling shifts the 
maser frequency antisymmetrically with respect to the detuning of the applied field about the Zeeman resonance . 

To second order in the Rabi frequency of the applied Zeeman field, \Xi2\ the small static-field limit of the maser 
frequency shift from the unperturbed frequency is given by 

Am- \Y \2(n0 .,0 N '^(7172 + l^24p)(7z/7b) ,„o^ 



\Xr. 



|2 



f uJc - ^2i\ Qclzi'^ + K) 



where 72 is the Zeeman decoherence rate, S = ujz — '-^23 is the detuning of the applied field from the atomic Zeeman 
frequency, K = ||-'i^24p/(7f + S'^), and pf^ — pgg ~ 7b/(27i) is the steady state population difference between states 
|1) and |3) in the absence of the applied Zeeman field. The first term in Eqn. ^results from the coherent two-photon 
mixing of the F = 1 levels as described above , while the second term is a modified cavity pulling term that results 
from the reduced linc-Q in the presence of the applied Zeeman field. We compared Eqn. ^ to experimental data 
from P-8, inserting the independently measured values of 1^241; 7fc7 7ii ^^'^ 72- By matching the fit to the data we 



extracted the Zeeman field parameters 



X12I and jz shown in Fig. ||. 
In addition to the shift given by Eqn. 23[ there is a small symmetric frequency shift due to the slight non-degeneracy 
of the two F ^ I, Amp = ±1 Zeeman frequencies. This term offsets the zero crossing of the maser shift resonance 
away from the average Zeeman frequency ^ (1^12 + 1^23), however the contribution is negligible at small static fields. 
Also, a reanalysis of the double resonance maser shift which included the effects of spin-exchange collisions 
showed that there is an additional hydrogen density-dependent offset of the zero crossing of the maser shift resonance 
from the average Zeeman frequency. Using the full spin-exchange corrected formula for the maser frequency shift , 
we calculated this offset and found that for typical hydrogen maser densities (n « 3 xlO^^ m~^), the offset varied 
with average maser power as approximately -50 /iHz/fW (assuming a linear relation between maser power and atomic 
density of ^ w 3x\'o°/m-3 )- described below, our masers typically have sidereal power fluctuations less than 1 
fW, making this effect negligible. 

The applied Zeeman field also acts to diminish the maser power, as shown in Fig. ^, and to decrease the maser's 
line-Q. By driving the F = 1, Amp — ±1 Zeeman transitions, the applied field depletes the population of the upper 
masing state |2), thereby diminishing the number of atoms undergoing the maser transition and reducing the maser 
power. Also, by decreasing the lifetime of atoms in state |2), the line-Q is reduced. A very weak Zeeman field of 
about 35 nG (as was used in our Lorentz symmetry test) decreases the maser power by less than 2% on resonance and 
reduces the line-Q by 2% (as calculated using Eqn. 6 of |^^). The standard method of determining the average static 
magnetic field strength is to scan the Zeeman resonance with a large applied field and record the power diminishment 
(such as that shown in Fig. ^, open circles). From the applied field frequency at the center of the power resonance, 
which typically has a width of about 1 Hz, the magnetic field can be found with a resolution of about 1 pG. 



8 



N 

X 



Si 



10- 



1X2401 : 

71 = 

12-- 

Iz-- 
Yb 



0- ^ 



-10- 



: 2.82 rad/s 
: 1.02 rad/s 
4.13 rad/s 
2.80 rad/s 
: 2.30 rad/s 
= 0.86 rad/s 



GO 



I I 




- 10 



- 5 



I 



--5 



--10 



I I I I I I I 

-3-2-10 1 2 3 

Zeeman detuning [Hz] 

FIG. 3. Double resonance maser frequency shifts. The large open circles (maser P-8) are compared with Eqn. ^ (full curve) 
using the parameter values shown. The values of |Xi2 1 and 7^ were chosen to fit the data, while the remaining parameters 



were independently measured as described in subsection IIIB. The experimental error of each measurement (about 40 /iHz) is 
smaller than the circle marking it. The solid square data points are data from the Lorentz symmetry test (maser P-28). The 
large variation of maser frequency with Zeeman detuning near resonance, along with the excellent maser frequency stability, 
allows the Zeeman frequency (~ 800 Hz) to be determined to 3 mHz in a single resonance (requiring 18 minutes of data 
acquisition). The inversion of the shift between the two is due to the fact that for the P-8 data (open circles), the maser 
operated with an input flux of |2) and |3) atoms, while for the P-28 data (solid points), the typical input of |1) and |2) atoms 
was used. Changing between these two input flux modes is done by inverting the dir ection of the static solenoid field, while 



maintaining a fixed quantization axis for the state selecting hexapole magnet (see Sec. IV D) 
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FIG. 4. Double resonance maser power diminishment. The open circles, taken with an applied Zeeman field strength of 
about 560 nG, represent typical data used to determine the value of the static magnetic field in the maser bulb. The filled 
circles are maser power curves with an applied field strength of about 80 nG. Our Lorentz symmetry test data were taken with 
a field strength of about 35 nG, where the power diminishment is less than 2%. 



IV. EXPERIMENTAL PROCEDURE 
A. Zeeman frequency measurement 



To measure the F = \, IS.mp = ±1 Zeeman frequency, we applied an oscillating field of about 35 nG near the 
Zeeman frequency. This field shifted the maser frequency by a few mHz (at the extrema) , a fractional shift of about 2 
parts per trillion. Because of the excellent fractional maser stability (2 parts in 10^^ over our averaging times of 10 s), 
the shift was easily resolved (see the solid data in Fig. ||) . As the frequency of the applied field was stepped through the 
Zeeman resonance, the maser frequency (of perturbed maser P-28) was compared to a second, unperturbed hydrogen 
maser frequency (P-13). The two maser signals at w 1420 MHz were phase locked to independent voltage controlled 
crystal oscillator receivers. The exact value of the receivers' outputs were set by tunable synthesizers, which were set 
such that there was a 1.2 Hz offset between them. The two receiver outputs were combined in a heterodyne mixer 
and the resulting 0.8 s period beat note was averaged for 10 s (about 12 periods) with a Hewlett-Packard Model HP 
5334B frequency counter. The full double resonance spectrum consisted of 100 such points. For each spectrum, 80% 
of the points were taken over the middle 40% of the scan range, where the frequency shift varies the most. 

Once an entire spectrum of beat period vs applied Zeeman frequency was obtained, it was fit to the function 

° Ai(1 + k)2+,52(1_^)2 Ai(l + K)2 + (^ + r)2(l-K)2 ' 

A5S , Ae{l + K) 



(Ai - ,52 + A4)2 + 452^1 Ai(1+k)2+,52(1_^)2 

to determine the Zeeman frequency. Here S ^ ly — vz the Zeeman detuning of the applied field away from the 
Zeeman frequency vz, n = ^4/(^1 + S^) is the analog of the parameter K from Eqn. and r = (1.403 x 10^^) x j/| 
is the small difference between the two Zeeman frequencies 1^12 and 1^23. The first term Aq is the constant offset 
representing the unperturbed beat period between the two masers. The second and third terms comprise the first- 
order symmetric maser shift (not included in Eqn. |2^ but described in the text above); these two terms nearly cancel 



at low static field where r vanishes. The final two terms account for the two shifts given in Eqn. 23 
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FIG. 5. Results from a Monte Carlo analysis. The horizontal axis represents the shift of the Zeeman frequency as determined 
by our fits of over 100 synthetic data sets, the vertical axis is the number within each shift bin. The width of the Gaussian fit 
to the data is 2.7 mHz, representing the resolution of a single Zeeman frequency measurement. 



For our spectra in maser P-28 with small applied field amplitude (solid square data points of Fig. ^), typical fit 
parameters were: Aq = 0.84550 ± 0.00001, Ai = 0.141 ± 0.005, vz = 857.063 ± 0.003, A3 = 0.006 ± 0.010, A4 = 
0.029 ± 0.003, A5 = (3.2 ± 0.1) xlO"'^, and Aq = (-1 ± 5) xlO'^. The uncertainty in the Zeeman frequency was 3 
mHz. Also, A3 and Aq, the amplitude coefficients of the residual first order effect and the cavity pulling term, were 
consistent with zero. 

With our known value of 7b = 0.86 rad/s, and our measured value of 72 — 2.77 rad/s (from line-Q), the above set 
of fit parameters were consistent with the reasonable values X24 = 2.14 rad/s, 7z = 2.36 rad/s, 71 = 2.88 rad/s, and 
X12 = 0.40 rad/s (since A3 had such a large error bar, the value of X12 was chosen such the ratio of the maser shift 
amplitude in P-28 to P-8, shown in Fig. ^ is equal to the ratio of the squares of X12 for P-28 to P-8). 

To determine the number of points and length of averaging that optimized the Zeeman frequency resolution, 
we recorded several spectra with 50, 100, and 150 points at 5 s and 10 s averaging. We also varied the "density 
distribution" of points, including spectra where the middle 40% of the scan contained 80% of the points and those 
where the middle 30% contained 80% of the points (thus increasing the number of points in the region where the 
antisymmetric shift varies the most). With each of these spectra, we ran the following Monte Carlo analysis 
after fitting each scan to Eqn. |2j, we constructed 100 synthetic data sets by adding Gaussian noise to the fit, with 
noise amplitude determined by the unperturbed maser frequency resolution of about 40 /iHz. Each of these synthetic 
data sets was fit and a histogram of the fitted Zeeman frequencies was constructed. The resolution of each spectrum 
was taken as the width of the Gaussian curve that fit the histogram (see Fig. |^). As the total length of the scans 
increased, the resolution improved and converged to a limit of around 2.5 mHz. While the resolution improved slowly 
with increased acquisition time, it would have eventually begun to degrade due to long term drifting of the Zeeman 
frequency. (As will be described below, we found that the Zeeman frequency exhibited slow drifts of about 10-100 
mHz/day). We therefore chose a scan of 100 points at 10 s averaging, for a total length of about 18 minutes for 
our Lorentz symmetry test spectra. The results from the Monte Carlo analysis for one of these spectra indicated a 
Zeeman frequency resolution is 2.7 niHz (see Fig. 0). 



B. Data analysis 

Our net result combines data from three runs. During each data run, the 18 minute Zeeman frequency scans were 
automated and run consecutively. After every 10 scans, 20 minutes of "unperturbed" maser frequency stability data 
was taken to track the maser's stability. Each run contained about 10 continuous days worth of data, and each set 
contained more than 500 Zeeman frequency measurements, taken at w 18 minute intervals. 

For each run, the long term Zeeman frequency data was fit to a function of the form 

fit — {piecewise continuous linear function) + Siyz,aCOs{u!sidt) + 6i>z.[3S^'n{LUsidt) (25) 
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where 5vz,a and 5vz,i3 represent the cosine and sine components of the sidereal sinusoid. The time origin of the 
sinusoids for all three runs was taken as midnight (00:00) of November 19, 1999. The subscripts a and (3 refer to 
two non-rotating orthogonal axes perpendicular to the rotation axis of the earth. The total sidereal amplitude was 
determined by adding Svz.a and Si'z.f3 in quadrature. During each run, the Zeeman frequency drifted hundreds of 
mHz over tens of days. The piecewise continuous linear function, consisting of segments one sidereal day in length, 
was included to account for these long term Zeeman frequency drifts. This function was continuous at each break, 
while the derivative was discontinuous. 

The result of this analysis, where the fitting function (Eqn. ^) was applied to the full data set, was found to be in 
good agreement with a second analysis, where each individual day of data was fit to a line plus the sidereal sinusoid 
and the cosine and sine amplitudes of each day were averaged separately and then combined in quadrature to find 
the total sidereal amplitude. 



C. Run 1 

The cumulative data from the first run (November 1999) are shown in Fig. ^(a) and the residuals from the complete 
fit (Eqn. 25) are shown in Fig. ^(b). The data set consisted of 11 full days of data and had an overall drift of about 
250 mHz. 

To avoid a biased choice of fitting, we allowed the location of the slope discontinuities in the piecewise continuous 
linear function to shift throughout a sidereal day. We made eight separate fits, each with the location of the slope 
discontinuities shifted by three sidereal hours. The total sidereal amplitude and reduced chi square for each is shown 
in Fig. 0. We chose our result from the fit with minimum reduced chi square. 
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FIG. 7. (a) Total sidereal amplitudes for the first run. The different points are from different choices of slope discontinuity 
locations, (b) Corresponding reduced chi square parameters. The minimum value occurs with a slope break origin of midnight 
(00:00) of November 19, 1999. 
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Run 


Si^z.a [mHz] 


(Jq [mHz] 


Si'z,i3 [mHz] 


a/s [mHz] 


1 


0.43 


0.36 


-0.21 


0.36 


2 


-2.02 


1.27 


-2.75 


1.41 


3 


4.30 


1.86 


1.70 


1.94 



TABLE II. Sidereal amplitudes from all runs. 



As noted above, the error bar on a single Zeeman frequency determination was about 3 mHz. However, when 
analyzing a smooth region of long term Zeeman data (about 1 day) we calculate a standard deviation of about 5 mHz. 
We believe this error bar is due mainly to residual thermal fluctuations (see Fig. p^ . 

For our choice of slope discontinuity with minimum reduced chi square [^, the cosine amplitude was 0.43 mHz ± 
0.36 mHz, and the sine amplitude was -0.21 mHz ± 0.36 mHz. The total sidereal amplitude was therefore 0.48 niHz 
± 0.36 mHz. 



D. Field-inverted runs 2 and 3 



In runs 2 and 3, the static solenoid field orientation was opposite that of the initial run to further study the double 
resonance technique and any potential systematics associated with the solenoid field. With the static field inverted, 
and therefore directed opposite the quantization axis in the state selecting hexapole magnet, the input flux consists 
of atoms in states |2) and |3) (rather than the states |1) and |2)). Thus, reversing the field inverts the steady state 
population difference {pii — p^g) of Eqn. |2^ and acts to invert the antisymmetric double resonance maser frequency 

shift m. 

Operating the maser in the field reversed mode degrades the maser performance and subsequently the Zeeman 
frequency data. With opposed quantization fields inside the maser bulb and at the exit of the state selecting hexapole 
magnet, a narrow region of field inversion is created. Where the field passes through zero, Majorana transitions 
between the different mp sublevels of the F ^ 1 manifold can occur. This can alter the number of atoms in the 
upper maser state {F = 1, mp — 0, state |2)), which diminishes the overall maser amplitude and stability. In the 
field-inverted configuration, the maser amplitude was reduced by 30%, and both the maser frequency and Zeeman 
frequency were less stable. In addition, the field-inverted runs were each conducted soon after a number of rather 
invasive repairs were made to the maser | p8| . Thus, the quality of the latter two data sets was somewhat degraded 
from the first run (see Figs, ^(a) and ^(a)). The overall drift was larger (nearly 800 mHz over about 10 days), and 
the scatter in the data was increased, as can be seen from the residual plots from these runs (Figs, ^(b) and ^(b)) 
which have been plotted on the same scale as the residuals from the first run (Fig. ||(b)). 

The latter two runs were also less suitable for the piecewise continuous linear drift model used in the first run. In 
that case, the large slope changes were coincidentally separated by an integer number of sidereal days; in the last two 
runs, the larger and more frequent changes in slope were not. Therefore, only certain selected sections could be fit to 
the same model (Eqn. ^5|), significantly truncating the data sets. Due to all of these factors, the sidereal amplitudes 
and the associated error bars were up to an order of magnitude larger for the field-inverted runs than the first run. 
All values are shown together in Table |l[ 



E. Combined result 



The final sidereal bound, combining all three runs, was calculated using the data in Table ||. First, the weighted 
averages of the cosine and sine amplitudes, Si'z,a and Svz.fS, were found using the standard formula for weighted 
mean 

/^(.fi)/(4), M 

and their uncertainties were given by 
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FIG. 8. (a) Run 2 data (December 1999), with solenoid current fluctuations subtracted. To the measured Zeeman frequencies, 
we added 894.942 Hz. (Note the sign reversal from run 1 to account for the inverted field), (b) residuals after fitting the data 
to Eqn. UB. 
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FIG. 9. (a) Run 3 data (March 2000), with solenoid current fluctuations subtracted. To the measured Zeeman frequencies, 
we added 849.674 Hz. (Note the sign reversal from run 1 to account for the inverted field), (b) residuals after fitting the data 
to Eqn. p3. 
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(27) 



The sign reversal due to the field inversion was accounted for in the raw data, before the data were fit. Thus, the runs 
are combined using conventional (i.e., additive) averaging. The final sidereal amplitude A was calculated by adding 

/ — 2 - 2 

the mean cosine and sine amplitudes in quadrature, A — y ^ + Sv^ p. We measure a sidereal variation of the 
F — 1, Amp — ±1 Zeeman frequency of hydrogen oi A — 0.49 ± 0.34 mHz. 

We note that since we are measuring an amplitude, and therefore a strictly positive quantity, this result is consistent 
with no sidereal variation at the 1-sigma level: in the case where Svz^a Svz,/3 have zero mean value and the same 
variance a, the probability distribution for A takes the form P{A) = Au"^ exp(— A^/2(T^), which has the most 
probable value occurring at A = a. 



V. ERROR ANALYSIS 



In addition to our automated acquisition of Zeeman frequency data, we continuosly monitored the maser's external 
environment. At every ten second step, in addition to applied frequency and maser beat period, we recorded room 
temperature, maser cabinet temperature, solenoid current, maser power, ambient magnetic field, and active Helmholtz 
coil current (see Sec. V A). 



A. Magnetic systematics 

The F = 1, mp = ±1 Zeeman frequency depends to first-order on the z-component of the magnetic field in the 
storage bulb. Thus, all external field fluctuations must be sufficiently screened to enable a sensitivity to shifts from 
Lorentz and CPT symmetry violation. The maser cavity and bulb are therefore surrounded by a set of four nested 
magnetic shields that reduce the ambient field by a factor of about 32,000. We measure unshielded fluctuations in 
the ambient field of about 3 mG (peak-peak) during the day, and even when shielded, these add significant noise to 
a single Zeeman scan, as illustrated in Fig. ^(a). Furthermore, the amplitude of the field fluctuations is significantly 
reduced late at night, which could generate a diurnal systematic effect in our data. 

To reduce the effect of fluctuations in the ambient magnetic field, we installed an active feedback system (see 
Fig. [Tc| ) consisting of two pairs of large Helmholtz coils (2.4 m diameter). The first pair of coils (50 turns) produced 
a uniform field that cancelled most of the z-component of the ambient field, leaving a residual field of around 5 mG. 
A magnetometer probe that sensed the residual ambient field was placed partially inside the maser's magnetic shields 
near the maser cavity. This probe had a sensitivity of s = 1.7 mG/V. Due to its location partially inside the magnetic 
shields, the probe was screened by a factor of about six from external fields, reducing the sensitivity to s' = 0.3 mG/V, 
and producing a differential screening of 5300 between the magnetometer probe and the atoms. The magnetometer 
output was passed into a PID servo (Linear Research model LR-130), which contained a proportional stage (gain G 
— 33), an integral stage (time constant T^ = 0.1 s) and a derivative stage (time constant = 0.01 s). The correction 
voltage was applied to the second pair of Helmholtz coils (3 turns) which produced a uniform field {p = 14 mG/V) 
along the z-axis to nullify the residual field and actively counter any field fluctuations. Neglecting the small effect of 
the derivative stage, the overall time constant of this system was given by r = ri(l -I- s' /pG) « 0.1 s, about 100 times 
shorter than the averaging time of our maser frequency shift measurements (10 s). 

With this system we were able to further reduce ambient field fluctuations at the magnetometer by a factor of 3,000. 
The resulting unshielded fluctuations were less than 1 /xG peak-peak. The field recorded by the partially screened 
magnetometer probe is shown in Fig. |l^. The noise on a single Zeeman scan was reduced below our Zeeman frequency 
resolution, as shown in Fig. pT|(b). During our Lorentz symmetry test, we monitored the field at the magnetometer 
probe and placed a bound of ~ 5 nG on the sidereal component of the variation. This corresponds to a shift of 
less than 0.2 ^Hz on the hydrogen Zeeman frequency, three orders of magnitude smaller than the sidereal Zeeman 
frequency bound measured. 

The magnetometer |3^] used in the feedback loop was a fluxgate magnetometer probe (RFL industries Model 101) 
which consisted of two parallel high-permeability magnetic cores each surrounded by an excitation coil (the excitation 
coils were wound in the opposite sense of each other). A separate pickup coil was wound around the pair of cores. 
An AC current (about 2.5 kHz) in the excitation coils drove the cores into saturation, and, in the presence of any 
slowly varying external magnetic field oriented along the magnetic cores' axes, an EMF was generated in the pickup 
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FIG. 10. Schematic of the active Helmhohz control loop. A large set of Helmholtz coils (50 turns) cancelled all but a residual 
^ 5 mG of the z-component of the ambient field. This residual field, detected with a fiuxgate magnetometer probe, was actively 
cancelled by a servoloop and a second pair of Helmholtz coils (3 turns). The servoloop consisted of a proportional stage (gain 
— 33), and integral stage (time constant = 0.1 s) and a derivative stage (time constant — 0.01, not shown). The overall time 
constant of the loop was about r = 0.1 s. 
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FIG. 11. (a) Zeeman scan without the active Helmholtz feedback loop. The noise on the data is due to the left and right 
shifting of the antisymmetric resonance as the Zeeman frequency shifts due to 3 mG ambient field fluctuations, (b) Zeeman 
scan with active Helmholtz control. Ambient field fiuctuations were reduced to less than 1 /uG. 
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FIG. 12. Residual ambient magnetic field, after cancellation by the active Helmholtz control loop, sensed at tfie magnetometer 
probe. Each point is a 10 s average. These three days worth of data depict a Sunday, Monday and Tuesday, with the time origin 
corresponding to 00:00 Sunday. From these data it can be seen that for three hours every night the magnetic noise dies out 
dramatically, and that the noise level is significantly lower on weekends than weekdays. Nevertheless, with the active feedback 
system even the largest fluctuations (1 /iG peak-peak) causes changes in the Zeeman frequency well below our sensitivity 
(AS = IfiG ^ Auz = mfiHz). 



coil at the second and higher harmonics of the excitation frequency. The magnitude of the time-averaged EMF was 
proportional to the external field. The probe had a sensitivity of approximately 1 nG. 

Any Lorentz violating spin-orientation dependence of the energy of the electrons in the magnetic cores would induce 
a sidereal variation in the cores' magnetization and could generate, or mask, a sidereal variation in the hydrogen 
Zeeman frequency through the feedback circuit. However, based on the latest bound on electron Lorentz violation 
psf (10~^^ GeV), the Lorentz violating shift would be less than 10""'^^ G, far below the level of residual ambient field 
fiuctuations. Also, the additional shielding factor of 5300 between the probe and the atoms further reduced the effect 
of any Lorentz violating shift in the probe electrons' energies. 

With the ambient field kept nearly constant near zero, the Zeeman frequency was set by the magnetic field generated 
by the solenoid, and hence by the solenoid current. We monitored solenoid current fluctuations by measuring the 
voltage across the current-setting 5 kf2 resistor with a 5 1/2 digit multimeter (Fluke model 8840 A/ AF). By measuring 
the Zeeman frequency shift caused by large current changes, we found a dependence of around 10 mHz/nA. When 
acquiring Lorentz symmetry test data, we measured long term drifts in the current of about 5 nA (see Fig. p^ , 
significant enough to produce detectable shifts in the Zeeman frequency. Thus, we subtracted these directly from 
the Zeeman data. We measured a sidereal variation of 25 ± 10 pA on the solenoid current, corresponding to a 
sidereal variation of 0.16 ± 0.08 mHz on the Zeeman frequency correc tion. This systematic uncertainty in the Zeeman 



frequency was included in the net error analysis, as described in Sec. V C 



B. Other systematics 

The maser resided in a closed, temperature stabilized room where the temperature oscillated with a peak-peak 
amplitude of slightly less than 0.5 °C with a period of around 15 minutes. The maser was contained in an insulated 
and thermally controlled cabinet, which provided a factor of five to ten shielding from the room, and reduced the 
fiuctuations to less than 0.1 °C peak-peak, as shown in Fig. |lj. By making large changes in the maser cabinet 
temperature and measuring the effect on the Zeeman frequency, we found a temperature coefficient of about 200 
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FIG. 13. Solenoid current during the first data run. Each point is an average over one fuU Zeeman frequency measurement 
(18 mins). Since the Zeeman frequency is directly proportional to the solenoid current, we subtracted these solenoid current 
drifts directly from the raw Zeeman data, using a measured calibration. We find a sidereal component of 25 ± 10 pA to that 
correction, corresponding to a signal of 0.16 ± 0.08 mHz on the Zeeman frequency. This systematic uncertainty has been 
included in our overall error analysis. 



mHz/ °C. We believe this frequency shift was due mainly to the resistors which set the solenoid current, which had 
100 ppm/ °C temperature coefficients. We monitored the cabinet temperature and placed a bound on the sidereal 
component of the temperature fluctuations at 0.5 mK, which would produce a systematic sidereal variation of 100 
/iHz on the Zeeman frequency, about a factor of 3 smaller than the measured limit on sidereal variation in Zeeman 

frequency. 

As mentioned in Sec. [II C , spin-exchange effects induce a small offset of the Zeeman frequency given by Eqn. 
from the actual Zeeman frequency [^ . This would imply that fluctuations in the input atomic flux (and therefore 
the maser power) could cause fluctuations in the Zeeman frequency measurement. We measured a limit on the shift 
of the Zeeman frequency due to large changes in average maser power at less than 0.8 mHz/fW. (Expected shifts from 
spin-exchange are ten times smaller than this level (Sec. QIC). We believe the measured limit is related to heating 
of the maser as the flux is increased). During long-term operation, the average maser power drifted approximately 1 
fW /day (see Fig. |l^). The sidereal component of the variations of the maser power were less than 0.05 fW, implying 
a variation in the Zeeman frequency of less than 40 /iHz, an order of magnitude smaller than our experimental bound 
for sidereal Zeeman frequency variation. 



C. Final result 



We measured systematic errors in sidereal Zeeman frequency variation (as described in Sees. VA and VB) due 
to ambient magnetic held (0.2 /iHz), solenoid field (80 /iHz), maser cabinet temperature (100 /iHz), and hydrogen 
density induced spin-exchange shifts (40 /iHz). Combining these errors in quadrature with the 0.34 mHz statistical 
uncertainty in Zeeman frequency variation, we find a sidereal variation of the F — 1, Amp — ±1 Zeeman frequency 
in hydrogen of 0.44 ± 0.37 mHz at the 1-ct level. This 0.37 mHz bound corresponds to 1.5 x 10~^^ GeV in energy 
units. 
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FIG. 14. Temperature data during the first run. Each point is a 10 second average. The top trace shows the characteristic 
0.5 °C peak-peak, 15 minute period oscillation of the room temperature. The bottom trace shows the screened oscillations 
inside the maser cabinet. The cabinet is insulated and temperature controlled with a blown air system. In addition, the 
innermost regions of the maser, including the microwave cavity, are further insulated from the maser cabinet air temperature, 
and independently temperature controlled. The residual temperature variation of the maser cabinet air had a sidereal variation 
of 0.5 mK, resulting in an additional systematic uncertainty of 0.1 mHz on the Zeeman frequency. This value is included in 
the net error analysis. 
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FIG. 15. Average maser power during the first data run. Each point is an average over one full Zeeman frequency mea- 
surement (18 mins). We measure a sidereal variation in this power at less than 0.05 fW, leading to an additional systematic 
uncertainty in the Zeeman frequency of 0.04 mHz, which is included in the net error analysis. 
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FIG. 16. Coordinate systems used. The (X,Y,Z) set refers to a fixed reference frame, and the (x,y,z) set refers to the 
laboratory frame. The lab frame is tilted from the fixed Z-axis by our co-latitude, and it rotates about Z as the earth rotates. 
The a and /3 axes, described in Sec. span a plane parallel to the X-Y plane. 



VI. DISCUSSION 



A. Transformation to fixed frame 



Our experimental bound of 0.37 mHz on sidereal variation of the hydrogen Zeeman frequency may be interpreted in 
terms of Eqn. ^ as a bound on vector and tensor components of the standard model extension. To make meaningful 
comparisons to other experiments, we transform our result into a fixed reference frame. Following the construction in 
reference |Q, we label the fixed frame with coordinates (X,Y,Z) and the laboratory frame with coordinates (x,y,z), 
as shown in Fig. ^ We select the earth's rotation axis as the fixed Z axis, (declination — 90 degrees). We then define 
fixed X as declination = right ascension = degrees, and fixed Y as declination ~ degrees, right ascension = 90 
degrees. With this convention, the X and Y axes lie in the plane of the earth's equator. Note that the a, (5 axes of 



Sec. [VE, also in the earth's equatorial plane, are rotated about the earth's rotation axis from the X,Y axes by an 
angle equivalent to the right ascension of 71° T longitude at 00:00 of November 19, 1999. 

For our experiment, the quantization axis (which we denote z) was vertical in the lab frame, making an angle x ~ 
48 degrees relative to Z, accounted for by rotating the entire (x,y,z) system by x about Y. The lab frame (x,y,z) 
rotates about Z by an angle Oi, where is the frequency of the earth's (sidereal) rotation. 

These two coordinate systems are related through the transformation 
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TABLE HI. Electron, proton and neutron experimental bounds. 



Then, vectors transform as biab — T b fixed, while tensors transform as diab — T dfixedT ^■ 

As shown in equation (^, our signal depends on the following combination of terms (for both electron and proton): 

63 = 63 - mdao - Hi2- (29) 

Transforming these to the fixed frame, we see 

^3 = bz cos x + bx sin x cos fit + by sin x sm fit, 
d3Q = dzo cos X + dxo sin x cos fit + dyo sin x sin ^lt, (30) 
H12 = HxY cos X + Hyz sin x cos ^It + -ffzx sin x sin fit, 

so our observable is given by 

^^3 = [bz - mdzo - Hxy) cos x 

+ [by — mdyo — Hzx)s\'n.xs\n.U.t (31) 
+ [bx — Tndxa ~ Hyz) sin x cos fit. 

The first term on the right is a constant offset, not bounded by our experiment. The second and third terms each 
vary at the sidereal frequency. Combining Eqn. ^ (for both c^ and p) with Eqn. ||, we see 

\Ai.z\' - Wy - ruedl^, ~ H^zx) + K " m,£yo " ^Ix)]' ^ (32) 
+ Wx ~ ^ed\^ - + (h\ - m^d\^ - H^^f 

Inserting x = 48 degrees, we obtain the final result 



' {b\ + V {b\. + fef,) ^ = (3 ± 2) X 10-27 GeV. (33) 
Our 1-sigma bound on Lorentz and CPT violation of the proton and electron is therefore 2 x 10^^^ GeV. 



B. Comparison to previous experiments 



We compare our result with other recent tests of Lorentz and CPT symmetry in Table |II[ Although our bounds 
are numerically similar to the those from the ^^^Hg/^'^'^Cs experiment, the simplicity of the hydrogen atom allows us 
to place bounds directly on the electron and proton; uncertainties in nuclear structure models do not complicate the 
interpretation of our result. The recent limit set by the torsion pendulum experiment of Adelberger et. al. ||20|] on 
electron Lorentz and CPT violation casts our result as a clean bound on Lorentz and CPT violation of the proton. 



C. Future work 



To make a more sensitive measure of the sidereal variation of the Zeeman frequency in a hydrogen maser, it will 
be important to clearly identify and reduce the magnitude of the long term drifts of the Zeeman frequency. Possible 
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sources of these drifts are magnetic fields near the maser bulb caused by stray currents in heaters or power supplies 
in the inner regions of the maser. Also, the scatter of the Zeeman data points, believed to be due mainly to residual 
thermal fluctuations, should be reduced. Both of these objectives could be accomplished by carefully rebuilding a 
hydrogen maser, with better engineered power and temperature control systems. 
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